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SUMMARY 

The K^irnian-Pohlhausen integral method, as applied to conpressihle 
laminar boundary layers, was simplified by an analysis similar to the 
inconpressible Holstein-Bohlen method. Althou^ this sinplification is 
helpful for many calculations, it is of greatest value when conditions 
at the „ edge of the boundary layer are known from ejcperimental measure- 
ments. The analysis was conducted under the assumptions of a Prandtl 
number of 1, zero heat transfer, 6uad a linear viscosity- tenperature 
relation; velocity profiles are approximated by a fourth-degree 
polynomial. 

Results are presented so that velocity and t^nperature profiles, 
moinentum and displacement thicknesses, and wall shear stress can be 
calculated for flows over two-dimensional bodies with arbitrary free- 
stream velocity distributions. The resvilts are also applicable to 
flows over three-dimensional bodies with axial symmetiy through the use 
of Mangier *s transformation. 


INTRODUCTION 

Present-day theory of conpresslble laminar boundary layers permits 
the aerodynamicist to make satisfactory boundary-layer calculations for 
flows with zero pressure gradients. Exact solutions of the boundary- 
layer equations for flows with streamwise pressure gradients, however, 
exist only for several specific pressure distributions (reference l) . 
Solutions of the lamineir-boundary-layer equations fo?: flows with arbi- 
trary pressure gradients are of ever increasing inportance because it 
has been suggested that at sufficiently hi^ altitudes and speeds lami- 
rar boundary layers can exist at very hl^ Reynolds numbers (refer- 
ence 2), and because the pressiire gradients encountered over most wings 
and bodies eire not approximated by any of the aforementioned exact 
solutions. 
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The lack of exact solutions for flows with arhltreiry streamwlse 
pressure gradients necessitates the use of approximate solutions. The 
EarmAn-Pohlhausen method has long heen found to 3 rield satisfactory solu- 
tions for inconpressible flows^ provided, that the separation point is 
not approached. Howarth (reference l) has recently applied this method 
to cong)ressihle flows. 

The solution of the incompressihle or conpresslhle houndary-layer 
equations hy the Kacmdn-Pohlhaiisen method is conplex because it involves 
the numerical solution of a cumbersome differential equation, and 
because it involves second derivatives of the free-stream. velocity, 
which cannot be determined with any accuracy if the free-stream velocity 
is known only from esperimental measurements. For incon 5 >ressible 
fltiids, these difficulties have been overcome by the Holsteln-Bohlen 
sinpllflcation (reference 3), which eliminates the second derivatives 
of the free-stream velocity and sin^illfies the differential eqviatlon. 

In the present report of an analysis made at the MCA Lewis labo- 
latory, the Holstein-Bohlen sinplification is applied to the KarmAn- 
Pohl hausen method for congiresslble fluids. The momentum- integral equa- 
tion is derived using the transformation of HowEirth (reference l) 
together with a more general viscosity-teTTTperature relation than that 
used in reference 1. The greater part of the analysis presented fol- 
lows along the lines of reference 1, but the eiisphasis heiein is placed 
on simplifying the final, results of the analysis so that they can easily 
be applied to any practlceil problem. 

Although the present analysis is carried out for two-dimensional 
flows, it can be applied to three-dimensional flows with axial symmetry 
by using the transformation of Mangier (see appendix A). TMder this 
transformation, a three-dimensional body with a given pressure gradient 
is transformed to a related two-dimensional body. The results of the 
present study can be applied to this two-dimensional body. 


SYMBOLS 

The following symbols are used in this report; 
a speed of sound 

C proportionality factor used in equation (l) 

Cp specific heat at constant pressure 

f (Mq) function defined by equation (28) 
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fx(^) fuaction defined ty equation (28) 

fgCx) function defined "by equation (28) 

f(il) function defined "by equation (20) 

h increment of length 

i summation index 

K arbitrary constant 

k thermal conductivity 

M Mach number 

n transformation variable defined by equation (9) 

Er Prandtl number, IJLCp/k 

p static pressure 

r radial coordinate 

S Sutherland's constant (216° F) 

t static teng)erature 

u,v velocity in x- and y-dlrections, respectively 

x,y cartesian coordinates measured along body and peipendlcular to 
body, reepectlvely 

'll. function defined by equation (25) 

X ratio of specific heats 

A boundary- layer thickness in x,n coordinate system. 

A?*" function defined by equation (l9) 

8 boundary-layer thickness in x,y coordinate system 

5* boundary-layer displacement thickness 

n n/A 
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0 function defined ty equation (l8) 

0 ‘boxindary-layer momentum thickness 

X Pohlhausen parameter 

H absolute "vdacosity 

V kinematic viscosity 

p mass density 

T shear stress 

transformed stream fimction 
stream function 
Subscripts; 

0 conditions at stagnation point 

1 conditions at outer edge of boundary layer 

s standard condition 

Tf condition at solid boundary . 


ASSDMPTIOTIS 

The following assxmg)tlons are made in addition to the usual 
bounda3cy-layer assun^jtlons in order to make the problem amenable to 
solution; 

(1) The Prandtl number of the fluid is equal to 1. It has been 
shown in reference 4 that this assun^tion has no appreciable effect on 
flat-plate skin- friction calculations if the fl\ald is air (for M<5), 
but the assunption does limit the accuracy of heat- transfer calcula- 
tions and teng)eratiare distributions. 

( 2 ) The heat transfer at the solid boundary is equal to zero. 

This assuD 5 )tion should be reasonable for many practical applications. 
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( 3 ) lEhe viscosity and ten®erature are related linearly "by the fol- 
lowing expression: 



( 1 ) 


Chapman and Euhesin (reference 5) have shown that solutions of the 
houndary-layer equations based on equation (l) agree well with reality 
for flat-plate flows if the constant C is determined hy matching 
equation (l) with Sutherland’s relation at the solid boundary so that 


C 



(ts + S) 
(iv + s) 


( 2 ) 


This assunption should also be reasonable for the case of flows with 
strearawise pleasure gradients because the wall tenperature t^ is con- 
stant as a consequence of assunptions (l) and (2). 


AHAtrSIS 

Gongressible laTninar-boimdary-layer equations . - The momentum equa- 
tion describing the flow in the boundary layer is 

^ , ^U 1 dp , 1 ^ ^u'N r,\ 

“ s ^ - p di + p ^ I** 

The equation of continuity is 

(pu) + ^ (pv) = 0 (4) 

The flow of heat is described by the energy equation, which is 



A particular solution of the energy equation can be obtedned by multl- 
.plying equation (3) by pu and adding the result to equation (5), and 
by letting the Prandtl number equal 1. This solution is 



+ Cp 


t = 


constant 


( 6 ) 


and corre^pnds to the case where the lieat transfer at the solid bound- 
ary is equal to zero. 
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The equation of state can he expressed as 

^ = _P_ JL 
Ps Ps ^s 


( 7 ) 


The Bernoulli equation, ■which applies at the outer edge of "the 
hounda3ry layer. Is 


is 

dx 


duj^ 

- Pl^l ^ 


( 8 ) 


Ti^nsformatlon of Howarth . - In reference 1 Howarth Introduces a 
transformation which, when applied to equation (3), yields an equation 
similar to the Inconpresslhle momentum, equation. In the present report 
the transformation Is modified sl.lghtly to Include the proportionality 
factor C In equation (l). The transformation 'varlahles are 



where n alters the scale In the direction normal to the surface. The 
derlvatl-ves are expressed as follows: 


a ^ ^ ^ 

1 

^ /p B 

"Sy “ t "iSn ^ 


(10) 


Before transforming equation (3) It Is convenient to Introduce a stream 
function ilr(x,y) which satisfies equation (4). This function Is 
defined hy 



pv = - Ps 


( 11 ) 
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In terms of the transformed coordinate system, a function j^(x,n) can 
he defined as 


1 

.2 


lirom equations (lO), (n) , ami (l2). 


u = 


V = - 



fee 4. ^ 




'Ps 



( 12 ) 


(13) 


and, as shown in reference 1, from equations (l), (?), (8), (lO), 
piTid (l3), the mamentum equation hecomes 


The temperature term in this equation can he expressed in the form 


t , ^1^ - 

ti 2Cp ti 





from which 








( 15 ) 


( 16 ) 


“With the exception of the term within the bracket, equation (l6) is 
identical in form to the Incompressihle momentum equation. 


Momentum- integreil equation . - The momentum- Integral equation is 
obtained by Integrating each term of equation (l6) from the solid 
boundary to the outer edge of the boundary layer (see reference l): 



d0 ^■'^1 

dx -L dx 



( 17 ) 
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vhere 



(18) 


(19) 


The auxiliary functions 0 A* are similar to the momentum thick- 

ness 0 aufl the displacement thickness 8* In the physical coordinate 
system. 

Mollified TCarman-PohThausen solution . - In order to solve equa- 
tion (17} ~hy t.hfi K~4.- rmATi -Pn hThau fieTi method it is necessary to introduce 
the variahle t) = n/A and to assume that u/uj^ is a polynominal in i]. i 

As in the method of Pohlhausen, this polynomial is taken to he of the 
fou3rth degree. The coefficients of the polynomial are determined from 
the following boundary conditions; 


n = 0; 


u = 0 


n = A; 



® ^3 




u = ui 


^ _ 


Gkius the following form is obta±n,ed for the velocity profile: 


u 

HI 


= f (il) = 2 ti - 2ti^ + 1]^ + X [j 11 (1-Tl)^ 


( 20 ) 


where the parameter X is defined as 


, _ a 2 dui / 2 \ 

^ — J 


( 21 ) 
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Transformation, to the tj- coordinate and siibsequent Integration In egua- 
tlons (is) and (19 ) yield 

(^) 

and It Is found that 

Equations (22) ^ (E3), and (24) could he applied to equation (l7) 
directly and the result would he a differential equation for X (ref- 
erence l) . A slnpler solution can he obtained hy applying the substi- 
tution of Holstein-Bohlen: 


Z = 02/^ g 

or^ from equations (22) and (2l) , 


(25) 


Z — j y- _ - ^ % 




( 26 ) 


Eqxoation (l7) is multiplied hy ©/u^Vq^, and equations (25) and (26) are 
substituted in the resulting e3q>ression to yield 


X i + i - 




Ww 


Then, from, equations (22), (23), and (24), 

i ' tit’ll 


( 27 ) 
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■wiiere 



(28) 


and 


Z = 


fi(x) 


d^ 

dx 



r=i 

2 



(29) 


Equations (27) and (29) can now "be integrated simultaneously, and X 
can "be obtained as a function of x "by a method of numericaO. integra- 
tion given in detail in appendix B. The functions fp(X) and f 2 (x) 
are graphically presented in figure 1 and are tabulated in table 
f(Mj[_) is tabxiLated in table II. Boundary-layer Variables are still in 

the transformed system of coordinates, however, and a transformation to 
the physical (x,y) coordinates is required. 


Inverse transformation to physical coordinates . - The distance 
normal to the surface can be obtained from equation (9); 



or 
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But, from equation (l5). 


= 1 + ^ 


and, from equation (20), 


where 




■H ^ 


\2 2 


S3 = - ^ - 2X 
4: - 12 ^ 3 ® 

« _ 5X2 5X ^ 

^■5 ” 9 ’*' 3 ^ 

ag = - 3X + 4 

X2 5X . 

ay = --g-+-^-4 

X2X - 
^ " 3 ^ 

(The functions are tabulated In table I and presented graphically 

In fig. 2.) 

from which 
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where A is ottained from equation (2l). The houndaiy- layer thickness 
is deteimined hy letting tj equal 1: 



(- 0 . 00011 ^ 


0.00941 + 0.4175) 


) 


(33) 


The momentum thickness is defined as 



The wall shear stress tt-,,- Is defined as 



A 



'pC 

Ps 



(36) 
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The sei)axatiorL point is defined as the point -where the shear s-fcress 
•vanishes, or, from equation (36), -where \ = -12. It is enphasized, 
howe-yer, tha-t -fche Earman-Pohlhausen solution has been found -to be quite 
inaccrura-fce in the neighborhood of a separation point. 


SUMMARY OF METHOD 


The results of -this s-budy -will be summarized In the order they are 
needed for the solution of a particular problem. It is assumed that eiU 

( ^'^1 ^ 

free- stream conditions ^u^^, 1^, t-j_, p^^l known functions, of x. 

The standard conditions (vg, Pg, tg) are selec-ted as mean -values of -fche 

free-s-fcream conditions. Hexfc X is de-termlned as a function of x by 
sol-vlng equations (27) and (29) numerically: 


•where 


^ ^ Ri(^) -(%) + ^2C>^^ 

dx Uj^ •- J 

fl(?^) 


Z = 




" 3i5 (= 
f(Mi) = 


X 5X 
" 3 " 144 


•xn X 5X 
“ 3 “ l44 


f 
)(■■* 





Mi^ - 4 


(27) 


(29) 


(28) 


2 + (r-l) Mi^ 

A suggested procedure for sol-ving equation (27) is presented in appen- 
dix B. All functions of X and % are tabula-bed in tables I and II . 


Once 
be found: 


X is known as a function of x,- the folio-wing functions can 



A = 


( 21 ) 
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(22) 



(23) 

and 



^ P T^S 

[l + Mi^ (-0.0001X2 _ 0.0094X + 0.4175^ 

(33) 


ts ^ P 

(34) 

6* == + A Mi^ (-0.0001x2 _ 0.0094X + 0.4175^ (35) 



(36) 

and, for assimed valuss of 0^ t;-^ 1, 


u 

^1 

= 2ti - 2 ti3 + tj4 + X [| Ti (l-il)^J 

(20) 



(30) 

The velocity ten 5 )erature profiles can then he found in terms 

physical (x,y) coordinate system hy the following relation: 

of the 


L i=2 

(32) 
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where 


\2 2 




X2 

--- 2X 


®-4= “ 12 3 " ® 


= - 


^+^ + 4 

9 3 


^8 


5X^ 

12 


- 3X + 4 


- 4 

6 3 


X2 

36 


+ 1 


SPECIAL CASES 


Solution neax a stagoatlon point . - In the immediate neighborhood 
of a stagnation point, up is equal to zero and the Incompressible 
solution must apply. !Ehis solution is (reference 3): 


Xq = 7.052 


Zo = 


0.077 

dup/dx 




(37) 




= -0.0652 


d^p/dx2 

dup/dx J 


Flat-plate solution . - In the case of the flat plate, X = 0 and. 
the conditions at the outer edge of the boundary layer are constant and 
equal to the standard condition(s) . Equation (27) therefore becomes 
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( 


= 0‘'^70 
^1 


and 


Z = 


0.470 X 

Ui 


From equations (25) and (34)^ the momentum thiclmess is 


0 = 0.685 



(38) 


(39) 


The boundary-layer thickness is found from equations (22), (34), 
and ( 33 ) to he 


5 




(40) 


The displacement thickness, from eqixations (23) and (35), is 
5* = 1.75^/^2" + 1.392 Mi^ 

and, from equation (36), the vail shearing stress is found to he 

^ 1»"1 f . 1 ^ 

S Vvix 


(41) 


" 2.918 Vvqx 


(42) 


Velocity and tenperature profiles can he found as functions of t] hy 
letting X vanish in equations (20) and (30), respectively. The rela- 
tion of T] and y can then he obtained from equation (32). The error 
introduced in the flat-plate solution hy assuming that the velocity 
profile is a fourth-degree polynomial can he deduced hy comparing these 
results vlth the resiilts of reference 5; 
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+ 1.372 


(¥) 


~ 3.02 



(44) 

(45) 


Equations (43) and. (45) apply for any Erandtl number, -wiiereas equa- 
tion (44) applies for Pr = 1 only. 

It can therefore he seen that the results of the present analysis 
are quite accurate for the case of zero pressure gradient. Eo definite 
conclusions for the nonzero pressure gradient case can he made, however. 


Lewis Fli^t Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 8, 1951 
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APPEKDIX A 


OIRAHSFORMATION OF MAUGLER 

Hie 1 i3Tn-inm- -'h mTnflR'r y-1 PLyfiT equations for tlrree- dimensional flews 
with axial symmetry are 


+ = 
ox oy 


p dx P ^ V 


Bx 


(P^o^) 



0 


PCJ) 






<A1) 

(A2) 


(A3) 


where the har is used to differentiate the three-dimensional from the 
two-dimensional quantities^ and rQ(x) defines the radial coordinate of 

the hody in a meridional plane. The cont inui ty equation is satisfied 
hy the following stream function; 




^0 P ^ - f 


(A4) 


^0 P ^ - 1 

The following transformation varlahles were introduced hy Mangier 
(unavailahle report) : 


= k 2 


0 


ro^(x) dx 


y = K ro(x) y 


(A5) 


where K is an arhitrary constant. Thus 
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Tq ^ 

^ -Kr ^ 

■wbere the prime indicates differentiation with reject to x. 
It is further assumed that 


P(x,y) = p(x,y) 


t(x,y) = t(x,y) 


l^(x,y) = n(x,y) 
p(x,y) - p(x,y) 

= g t(x,y) 


Equations (A1), 
form: 

(A2), and (A3) can then he ti’ansformed to the 

following 



3u , Su 1 dp , 1 S 

(A6) 




^ (A7) 


PCp 

/ Bt 3t\ ^ ^ 

(A8) 

where 


u(x,y) = u{x,y) 

v(x,y) - |j(5,y) - u(5,y^ 

(A9) 


Eqtiations (A6)^ (A7), and (A8) are now identical to equations (3), (4), 
and (5). The solution of these equations, as presented in the text of 
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■bhis report, ‘ttierefore applies "to 1ihe •bhree-dlTnenslonal problem, "but 
the final results must he r^transformed to the three-dimensional quan- 
tities. The boundary- layer thickness .5 Is not defined In a solution 
of the differential boundary-layer equations, but Its definition Is 
required for the solution of the Integral equation. The transformation 
of 6 can best be accon 5 )llshed by examining the velocity profiles: 



J_ = 

^1 ^l\^J 


Tram, equation (A9), however. 


u u 



It therefoie follows that 


3 : 

5 6 


•vdience, from equation (A5), 


6 = ETq(x) 6 (AIO) 


The momentum thickness Is transformed as follows; 



= ^ ^Cx) 


(All) 
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and similarly, 

5*(x) 


where the quantities without the bar (representing two-dimensional 
quantitlesj are given by equations (33) to (36). 

Therefore, if it is desired to determine any of the boundary-layer 
characteristics at a distance x from, the nose of an axially symmetric 
body, the corresponding two-dimensional distance x is first found from 
equation (A5). The two-dimensional quantities (&, B*, 0, t^, and so 

forth) are then found using equations ^33) to (36), and are transformed 
to the three-dimensional quantities at x throu^ the use of equa- 
tions (AIO) to (a 12). The t 2 ^ansformation app3J.es -for flows with or 
without streamwise pressure gradients. 




0 


= ^0 


] 


(A12) 
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APPEKDIX B 


SOLOTIOK OP KEFFEBEKTIAL EQU/VELOK BI RUHGE-EIIETA MEOHOD 

The Runge-Kutta method of finite differences (reference 6) has 
teen found to he very satisfactory for solving equation (27 ) . It is 
believed that a detailed discussion of the method, as applied to the 
present problem, 'will be of help to the reader. 


It is sviggested that a table of the following form be used: 



The distance along the body (x) is subdivided into a number of incre- 
ments h, and the free- stream variables Uq_, dU 2 /dx, and at each 

point are tabulated in columns (2), (3), and (4). Column (5) follows 
from equation (28), and column (6) is easily tabulated. Colu m n s (l) 
to (e) can be tabulated for all x, whereas for columns (?) to (ll), 
each row must be calculated separately. 

The value of Zq = (column (7), row l) is given by eqiia- 

tion (37) for stagnation-point flow, or is equal to zero if ©q = 0. 
The value of fQ^(X) (colimm 8) is obtained from equation (29), or in 
column notation. 


fl(x) = (7) (3) (6) 
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ColiJinns (9) and (lO) are obtained from figure (l). Column (u) 
lists tbe Runge-Kutta parameter k, wMch in column notation is 

k = ((8) (5) + (10)) b 

Rows 1 to 4 can now be con^jleted in this manner. At the conrpletion of 
row 4 tbe total Increment in Z can be computed using tbe following" 
equation: 

AZ = I (ki + 2k2 + 2k3 + k^) 

Tbe entire process is tben repeated for rows 5 to 8. Tbe value of Z^ 
(column ( 7 ), row 5) is equal to Zq + AZ. Tbe procedure is carried out 

for a1 1 increments of and as a result X is obtained as a function 
of X. 
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TABLE I - FUNCTIONS OF X 


DO 




1 



»2 


^4 

^5 

^6 

®7 


- 0 . 0001X2 

- 0 . 0094X 

40.4175 

a*/a 

e/A 

2 + x/6 

12 

0.0948 

0.1422 




- 56.000 

28.000 

- 8.000 

1.000 

0.2888 

0.2000 

0.0889 

4.000 

11 

.0941 

.1426 

14.695 

- 42.167 

57.084 

- 44.889 

21.417 

- 5.833 

.694 

.3008 

.2083 

.0925 

3.833 

10 


.1438 

13.444 

- 36.667 


- 34.889 

15.667 

- 4.000 

.445 

.3125 

.2167 

.0958 

3.667 

9 

.0882 

.1461 



37,750 

- 26.000 

10.750 

- 2.500 

.250 

.3240 

.2250 

.0990 

3.500 

8 


.1495 

fflnti 

- 26.667 

29.334 

- 18.221 

6.667 

- 1.333 

.111 

.3353 

.2333 

.1019 


7 

.0767 

,1544 


- 22.167 

21.750 

- 11.556 

3.417 

-.500 

.028 

.3463 

.2417 

.1047 


6 

.0689 

.1608 





1.000 

0 

0 

.3572 

.2500 

.1071 

3.000 

5 

.0599 

.1687 

8.028 

- 14.167 

9.084 

- 1.555 

-.583 

.167 

.028 

.3678 

.2583 

.1094 

2.833 

4 

.0497 

.1784 

7.111 

- 10.667 


1.778 

- 1 . 333 - 

0 

.111 

.3782 

.2667 

.1115 

2.667 

3 

.0385 

.1898 

6.250 


-.250 


- 1.250 

-.500 

.250 

.5883 

.2750 

.1133 

2.500 

2 

.0264 

.2030 

5.445 

- 4.667 

- 3.667 

5.111 

-.333 

- 1.333 

.445 

.3983 

.2833 

.1149 

2.333 

1 

.0135 

.2181 

4.694 

- 2.167 

- 6.250 

5.111 

1.417 

- 2.500 

.694 

,4080 

.2917 

.1163 

2.167 

0 

0 

.2350 

4.000 

0 

- 8.000 


4.000 

- 4.000 

1.000 

.4175 

.3000 

.1175 

2.000 

-1 

-.0140 

.2536 
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1.853 
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-3 
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2.250 

.4447 

.3250 

.1197 

1.500 

-■4 

-.0575 


1 . 77a 

5.333 

- 6,667 

- 11.556 

22.667 

- 13.333 

2.778 

.4533 

.3333 

.1199 

■ 1.133 

-5 

-.0720 

.3450 

1.361 

5.833 

- 4.250 

- 18.221 

29.417 


3.361 

,4617 

.3417 

.1200 

1.167 

-6 

-.0862 

.3717 



- 1.000 

- 26.000 



4.000 

.4699 

.3500 

.1198 

1.000 

-7 

-.0999 

.3992 

.695 

5.833 

3.084 

- 34.889 

45.417 

- 23.833 

4.694 

.4778 

.3583 

.1195 

.833 

-8 


.4279 

.444 

5.333 

8.000 

- 44.889 

54.667 


5.445 

.4855 

.3667 

.1189 

.667 

-9 

-.1255 

.4576 

.250 


13.750 

- 56.000 



6.250 

.4930 

.3750 

.1181 

.500 

-10 

-.1369 

.4874 

.111 

3.333 

20.333 

- 68.222 

75.667 

- 37.333 

7.111 

.5003 

.3833 

.1170 

.333 

-11 

-.1474 

.5181 

.028 

1.833 

27.750 

- 81.556 

87.417 


8.028 

.5074 

.3917 

.1158 

.167 

-12 

-.1567 

.5486 

0 

0 

36.000 

- 96,000 



9.000 

,5142 

.4000 

.1143 

0 
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TABLE II - FUNCTIONS OF 


0 0 

.2 


1 a. W 2 

1 + _ Ml 

f(Ml) 

■ 1.000 

-2.000 

1.008 

-1.974 

1.032 

-1.861 

1.072 

-1.698 

1.128 

-1.489 

1. 200 

-1.250 

1.288 

-.994 

1.392 

-.733 

1.512 

-.476 

1.648 

-.231 

1.800 

0 




0.968 

1.152 

1.352 

1.568 

1.800 

2.048 

2.312 

2.592 

2.880 

3.200 


1.968 

2.152 

2.352 

2.568 

2.800 

3.048 

3.312 

3.592 

3.880 

4.200 


f(Mi) 
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Figure 2. - Concluded. Auxiliary functions uss(!^ In equations (51) and (52). 
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